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We demonstrate that at a filling of n = 1.5, an hexatic insulating state obtains in the extended 
Hubbard model on a triangular lattice. Composed of two tetragonal sublattices with fillings of n = 1 
and n — 2, the insulating state is charge ordered and possesses an antiferromagnetic superlattice 
with dimension a x v3- Two distinct energy scales arise in our model, a charge gap for the insulator 
and the effective exchange interaction in the antiferromagnet. Our model is capable of explaining 
both qualitatively and quantitatively the Hall coefficient including the sign change, the temperature 
dependence of the resistivity and the persistence of antiferromagnetism above the insulating state. 



The cobaltates [lj], Na x CoC>2 are layered anisotropic 
materials that share many of the complexities of the high- 
temperature copper-oxide superconductors. They con- 
sist of insulating layers of Na + ions separated by C0O2 
planes. As in the cuprates, superconductivity in the hy- 
drated cobaltates[l| , Na x Co02-yH 2 0,&T:ises from doping 
a half-filled Mott insulating state. In the cobaltates, the 
half-filled band arises in the parent material, NagCo02, 
from a complete filling of each of the e g (t2 g ) orbitals while 
only a single electron resides in the d z 2 state @, How- 
ever, the underlying lattice is triangular as opposed to 
the square lattice of the cuprates. The concentration, 
x, measures the number of Co" 1 " 4 ions converted to Co +3 
ions for every Na atom added. Superconductivity[H in 
Na x Co02 ■ yH^O, is confined to a dome-like region as 
a function of the effective cobalt valence, which, in this 
case, is determined^ by the net hydronium and sodium 
ion contents 

While the cobaltates are generally metallic, they ex- 
hibit a novel insulating state[5| at x = 0.5. At present, 
the insulating behaviour at x = 0.5 remains a mystery. 
We propose that the key mechanism driving the insulat- 
ing state is the elimination of magnetic frustration by the 
formation of a charge-ordered state. Two distinct energy 
scales exist in the charge-ordered state: 1) the effective 
in-plane antiferromagnetic exchange coupling, </, and 2) 
the energy gap in the insulator, A. We show that for 
reasonable system parameters, A < J and hence the in- 
sulating state obtains at a lower temperature than does 
the antiferromagnetic (AF) transition. Our model pre- 
dicts electron-like transport below the AF transition but 
hole-like transport above, as is seen experimentally [Bj]. 

There are now distinctly different probes of the insu- 
lating state at n — 1.5. Transport measurements^ find 
that the resistivity at n = 1.5 is temperature indepen- 
dent from 60 < T < 300K but undergoes a sharp upturn 
as the temperature decreases below 53AT. Further, the 
Hall coefficient changes sign from positive to negative 
as the temperature is lowered through 88K. Unlike the 
square lattice which is particle- hole symmetric at n = 1, 
such is not the case for the triangular lattice. Hence, the 
vanishing of the Hall coefficient in the Na^CdC^ requires 
a non-trivial electronic state. The magnetic susceptibility 



X exhibits a slight kink at 53if, implying a phase tran- 
sition. In addition, neutron diffraction (^,[3] and NMRQ 
experiments demonstrate that the Na atoms form an or- 
thorhombic superlattice with lattice vectors (a\/3x, 2ay). 
Sodium atom ordering has been promoted as the a cen- 
tral ingredient in driving the insulating state. Consistent 
with such ordering is an enhancement of the thermal 
conductivity || below 53K only at n = 1.5. However, 
recent experiments on hydrated _/Vao.4iCo02 hi which 
the conductivity has also been observed to diverge be- 
low 53K call this view into question. Electrochemical 
measurements [9j indicate that charge transfer between 
the oxonium ion (H^O + ) and Co results in a formal va- 
lence of Co of 3.5 rather than the nominal value of 3.59 
obtained from the sodium content. Consequently, the 
authors conclude that the key ingredient in driving the 
insulating state is the presence of Co with a formal va- 
lence of 3.5 not sodium ordering. 

To understand what is special about the filling of n — 
1.5 on a triangular lattice, consider the square lattice 
at n = 1. At n = 1, the Coulomb repulsion energy is 
minimized and all magnetic frustration is eliminated on 
a square lattice by forming an antiferromagnet. Such 
a state is insulating because of the charge gap, U, for 
doubly occupying the same site. Our key assumption 
for the triangular lattice is that minimizing the Coulomb 
energy and eliminating Co magnetic frustration drives 
the insulating state at n = 1.5. At a filling of n = 1.5, 
all Coulomb repulsion energy and magnetic frustration 
can be eliminated by filling the lattice with alternating 
rows of doubly and singly occupied sites. Consequently, 
we consider the state in Fig. ([T]) in which the doubly and 
singly occupied sites occupy the corners of orthorhombic 
unit cells. 

Such a state can be thought of as a triangular lat- 
tice split into two tetragonal lattices, one for the double 
occupancy and the other for single occupancy with anti- 
ferromagnetic correlations. In fact, in terms of a unit cell 
defined by the basis vectors, (1,0) and (l/2,V3/2), the 
state proposed here preserves translational symmetry but 
breaks rotational invariance by 2n/3. Hence, the state 
is a hexatic nematic phase similar to that proposed for 
x = 1 in the Emery model [Io| in which nearest neighbour 
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FIG. 1: Charge ordering of Nao.sCo02- The circles repre- 
sent cobalt atoms. Notice that the triangular lattice is bro- 
ken into two orthorhombic lattices one with the double oc- 
cupancies and one with antiferromagnetic ordering. Moving 
an electron from a doubly occupied site to a singly occupied 
site will increase its energy by J where J = 4t 2 /U is the 
antiferromagnetic coupling constant. A next-nearest neigh- 
bour superexchange interaction, Ji , also exists in the vertical 
direction. The diagonal hopping process which leaves antifer- 
romagnetism intact is also depicted. 



Coulomb interactions play a pivotal role. On a triangu- 
lar lattice with n = 1.5, there is another ordered ground 
state, the zig-zag state[ll[, which can be obtained by 
a zero-energy cost rearrangement of the cobalt ions in 
the charge-ordered state. Consequently, which state is 
favoured experimentally is determined by factors such as 
the lattice structure and the sodium ordering pattern. 
Nonetheless, to establish that either the charge-ordered 
or zig-zag states is an insulator necessitates an analysis 
of the distinct transport channels, which we now present. 

We analyze the energetics for transport in the context 
of the extended Hubbard model, 



H = -ty^c| CT Cj CT + h.c. + U n^riii 

(ij) i 



V 



where Ci a {c\ a ) annihilates (creates) an electron with spin 
a on site i, [/, the energy cost for doubly occupying the 
same site and V the nearest-neighbour Coulomb repul- 
sion. The importance of the nearest-neighbour interac- 
tion in stabilizing the insulating state can be illustrated 
as follows. Since along each singly occupied row, charge 
transport is blocked as a result of the charge gap, U, it 
suffices to show that Fig. ([TJ is an insulator by demon- 
strating that the doubly occupied sites are immobile. 
Any electron in a doubly occupied site is surrounded by 
an equal number of spin up and spin down electrons in 
the neighboring sites, resulting in a net zero interaction 
energy. However, if an electron hops to a singly occupied 
site, it will be surrounded by three double occupancies, 
two electrons of the opposite spin, and one of the same 
spin. The net increase in the energy will be V + J where 
J = 4f 2 /[/[12} ■ We will call such a process as the creation 
of a doublon (on a singly-occupied row) and holon (on a 
doubly-occupied row) pair which costs an energy V + J. 



Once, a doublon-holon pair has been created, there are 
four distinct transport mechanisms: 

1) Holon motion along the otherwise doubly occupied 
superlattice: Transport of the holon vertically along the 
superlattice requires two-step hopping which does not 
cost any energy in the intermediate state. The ampli- 
tude for such a process is t. 

2) Doublon motion in the singly occupied row: It costs 
an energy nJi, were n is the number of hops. Hence, a 
doublon is linearly confined and does not propagate along 
a chain. 

3) Transport of a doublon along a diagonal through 
a series of concerted two-step hoppings as illustrated in 
Fig. ([T]): Doublon hops along the diagonal, however, do 
not generate any further energy costs and hence consti- 
tute the only truly 2-dimensional transport in the an- 
tiferromagnetic phase. As doublon motion propagates 
an electron, Rh < 0. This type of diagonal hopping 
process which leaves the antiferromagnet intact requires 
both the triangular lattice and the charge ordered state 
in Fig. ©. 

4) Holon transport along either one of the diagonals. 
Such hopping depends on the surrounding spin. 

On the other hand, in the absence of antiferromag- 
netic correlations, both the holon and doublon hop in 
their corresponding superlattices with the same hopping 
matrix elements but with opposite sign. All such trans- 
port mechanisms cease once the holon and doublon are 
bound. The criterion for holon-doublon binding can be 
derived from a 2-dimensional model, 

H = E 4 ( d Uj + h.c.) + tij civ + h.c.) (2) 

(*,i> 

+ E ( y + J)( h l*hi° + 4di) ~tJ2(4hl + h.c.) 

i i 

which includes all of the distinct hopping processes delin- 
eated above. The operator d k creates a doublon, while h\ 
creates a spin holon on the singly-occupied row. tfj and 
are the hopping matrices for the doublon and holon, 
respectively: 
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where a = ((— l) ri Si) m for some positive m is a function 
of the ordering parameter of the singly-occupied super- 
lattice which is temperature dependent and J = aJ . In 
our calculation, a has the following form 
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and constitutes the only free parameter in our compari- 
son with experiment. Both a and the modified form of 
J reflect the fact that the hopping processes depend on 



the antiferromagnetic ordering. As a consequence, J en- 
ters the energy cost only if antiferromagnetism is present. 
A i:) = t(5 ix ± lt j x 5 iy>jy + 8 ixtj:c 5 iy ±i tjy ) corresponds to the 
hopping mechanism for both the doublon and holon in 
the absence of antiferromagnetism. If the electron on 
the singly-occupied superlattice is antiferromagnetically 
correlated, the hopping mechanism for the doublon and 
holon becomes 



(5) 
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Likewise, a enters the diagonal hopping term be- 
cause such processes only enter if antiferromagnetism is 
present. Noting that the creation of a doublon-holon pair 
on two adjacent rows costs a total energy 2(V + J), we 
can write the non-interacting energies of the doublon and 
holon as 



e d , h (k) = V + J + t d ' h {k), 



(6) 



with t d (k) — —2at(cos(k x + V3k y ) + cos(k x — ^/2>k y )) — 
2£(l-a)(cos(fc x ) + cos(V3fc y )) and t ha (k) = 2t(cos(k x ) + 
cos(i/3ky)) + atcos(k x + a^/Sky) are the Fourier trans- 
forms of tfj and t^, respectively. Localization corre- 
sponds to a binding of a spinon and a holon. The solu- 
tion to this model, obtained via a Bogoliubov rotation, 
has the quasiparticle spectrum, 

E~(k) = + cos 2 9 k e d {k) - sin 2 9 k e h (k) -ism 20* 
E+(k) = - sin 2 9 k e d (k) + cos 2 9 k e h (k) -tsm29 k (7) 

where tan29 k = — 2t/(sd(k)+£h(k)) which corresponding 
to negative and positive charge excitation respectively. 
As is clear from Fig. ([2]), this model admits a metal- 
insulator transition with a well-defined gap, A, between 
the localized and the extended states. This gap is zero 
for V + J — At < but otherwise given by 



Aw (V + J)(l- 
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V + J 



(8) 



Experimentally, J w lOmeV and t — 0.2eV[l^; conse- 
quently, there is no solution to Eq. © if V = as in the 
Hubbard model. However, V is not known. Hence, to de- 
termine A, we must rely on experiments. To describe the 
experimental upturn of the resistivity at 53K , A w 53-KT, 
which results in V ~ 880mey > At. This value of V is 
reasonable because even if we use the reduced value of 
the on-site energy, U/V3 = 4.0eV/1.73 = 2.3ey[l3 re- 
ported recently for the insulating state, we find that V is 
of order [7/3. The importance of this estimate is that V 
is large enough to stabilize the charged ordered insulat- 
ing phase. We computed the resistivity using the Kubo 
formula within the diagonalised quasiparticle spectrum 
in Eq. ([7]). Fig. ((2^) shows that the computed resistiv- 
ity adequately reproduces the experimental trends. . The 
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FIG. 2: (a) Resistivity (dashed line) and band structure for 
the model in Eq. J2J. Two bands corresponding to negative 
and positive charge quasiparticles with a gap A are shown, 
(b) Hall coefficient for Eq. ([2]). The experimental data for 
both the resistivity (squares) and Hall coefficient (circles) are 
taken from Foo, et. al.0]. 



only free parameter is m in a which we set to 4. We have 
found, however, that the resistivity is insensitive to m as 
the divergence is governed primarily by the gap which is 
set by experiment. A key consequence of the insulating 
state proposed here is the emergence of an antiferromag- 
netic superlattice with dimension a x y/3a on the singly 
occupied sites. In fact, this is our key experimental pre- 
diction. This can be established rigorously. There is a 
weak antiferromagnetic coupling, J\ between electrons 
separated by \[Za on neighbouring singly occupied rows 
which is mediated by intervening doubly occupied sites. 
For this exchange interaction to be well-defined, doublons 
should not be created; hence V + J > At. Consequently, 
J i should be a function of V . 

To determine the rough magnitude of the effective 
spin coupling constant, Ji, in the sublattice, we con- 
sider the A th order perturbation process which flips the 
spins with respect to the charge ordered state and ob- 
tain J± = A(t/V) 2 J. For the experiments at hand, 
V ~ .886^ w At, implying that J\ w J. Furthermore, 
the effective horizontal spin coupling constant, J, should 
be renormalilzed by the same A th order process such that 
J — > J' = J + Ji m 1.25 J. Through the identical mech- 
anism, next-nearest neighbour sites (that is, sites along 
the diagonal) in the singly occupied sublattice are also 
antiferromagnetically coupled. However, in this case, 
only one path exists. Consequently, the next-nearest- 
neighbour exchange interaction is J2 = Ji/2. Hence, the 
singly occupied sublattice obeys the Heisenberg Hamil- 
tonian 

-^singl.occup. — J ^ • Sj-i-^ T J\ S^Sj-)_y 

% i 
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where x(y) denotes a nearest-neighbour along the x-axis 
(y-axis) . In the absence of J 2 , the system develops long- 
range antiferromagnetic order for temperatures less than 
J off = \/JJi L3ro< \ ' 1 7 . For J 2 ^ 0, antiferromag- 
netism obtains as long as Ji < J c ff/2[16] which is satis- 
fied in our case. As the effective low-energy description 
of the insulating state has now been reduced to a half- 
filled rectangular lattice (by virtue of integrating out the 
doubly occupied sublattice), particle- hole symmetry has 
now been reinstated. 

Consequently, the Hall coefficient should vanish as 
T — > for n — 1.5 as is observed experimentally 
Once the transition to the AF phase occurs, both the 
non-interacting quasiparticles with negative and positive 
charge can contribute to the Hall coefficient as men- 
tioned previously. But, according to Fig. [2j the nega- 
tively charged quasiparticles have a lower energy. Hence, 
Rh < below the AF transition. Above Tjy, no or- 
der exists but the stripe is still stable. In this regime, 
the sign of the Hall coefficient is determined strictly 
by the filling which for n = 1.5 should be hole-like. 
Since the Hall coefficient vanishes at T = 0, it is pro- 
portional to the charge difference between particles and 
holes. Hence, we computed the Hall coefficient by in- 
tegrating 2(f(E + (k)) — f(E~(k))/e over all momenta 
where / is the Fermi function. As Fig. (|2| indicates, 
the Hall coefficient contains has a sign change and a 
minimum at T « 0.5J that are in superb quantitative 
agreement with experiment. Likewise, the resistivity has 
a plateau where Rjj acquires its minimum value, as ob- 
served experimentally Q). Here again we used m — 4 and 
found that the overall qualitative features are insenstivive 
to the choice of to. The robustness of the qualitative fea- 
tures of Rh and the resistivity lend credence to Eq. ((2]) 
as the reductive model for the insulating state. 

Within mean- field theory [l7|, the upper bound to the 
temperature at which there is a transition to this AF 
state is given by T c = 2J c g. For reasonable system pa- 
rameters, this energy scale exceeds that of the gap for 
the insulating state and hence is consistent with recent 
neutron scattering experiments [3] which find charge or- 
dering and antiferromagnetism consistent with Fig. ([1]) 
above the temperature at which the resistivity exhibits 
an upturn. The relative ordering of the energy scales is as 
follows: the energy scale for destroying the stripe order 
is set by the nearest-neighbour Coulomb energy V. An- 
tiferromagnetic order is destroyed above J e g. Between 
A and J c g, AF order exists but the holons and dou- 
blons are mobile. Antiferromagnetic order in the plane is 
maintained because diagonal doublon hops do not gener- 
ate ferromagnetic bonds but rather exchanges doublons 
between the two sublattices. This ultimately mixes the 
valence between the cob alt atoms on the singly and dou- 
bly occupied rows. Hence, rapid charge fluctuations on 
a time scale set by Ti/t ~ 10~ 14 s could be the origin of 
the absence of two distinct Co environments reported in 



recent NQR experiments [19j. Nonetheless, such mixing 
ceases in the insulating state as the holon and doublon are 
bound. Hence, our T = state contains two distinct Co 
environments. High field NMR Co NMR experiments [iol] 
find, contrary to the NQR experiments, two distinct Co 
environments corresponding with effective moments con- 
sistent with the charge ordered state proposed here. 

To summarize, in NCI0.5C0O2, the ordering^, [il| of 
the sodium ions induces charge ordering in the CoOi lay- 
ers so that both the double occupancies and the singly 
occupancies reside on different tetragonal lattices. How- 
ever, charge localization obtains not from pinning by the 
sodium superlattice but by a delicate balance between 
local Coulomb repulsions and the kinetic energy in the 
C0O2 plane. The key prediction of this work that AF 
order on a superlattice with dimension a x ^/3a has been 
confirmed experimentally [3]. Further, the diagonal hop- 
ping mechanism identified here illustrates that the sign 
change of the Hall coefficient at Tn and the apparent 
itineracvfl9| of the antiferromagnetism are intimately re- 
lated. 
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